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SIKULATIOfI ANALYSIS OF A MlCP.OCC“PUTER-EASED 
LOU-COST OMEGA NAVIGATION SYSTEM 
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ABSTRACT 


Previous reported research has resulted In the develop- 
ment of a very low-cost Omega sensor processor [1,2, 3, 4], 
Currently under investigation is the implementation of a 
low-cost navigation processor and its Interrace with the 
sensor to complete the hardware-based Ohio University Cf!S. 
The iundatrental concept under investigation in this study 
is: If navigation processor functions can be performed by 

an inexpensive microcomputer , how many of the sensor pro- 
cessor functions can also be handled by innovative soft- 
ware? 

To explore this concept, computer simulation of sensor 
processor functions is underway. An input data base of 
live Omega ground and flight test data has been created. 
The Cr.cga senscr anc microcomputer interface '■cdules used 
to collect the data are functionally describes. Automatic 
synch.ronizaticn to the Omega transmission pattern is ces- 
crlbed as an example of the algorithms developed using 
this data base. 


The current status of research on a proposed micro- 
computer-based, low-cost Cmeca flavigaticn System (C:;S) is 
described. The design approacn emphasizes minirum hard- 
ware, maximum software, and the use of a low-cost, cermer- 
cially-available nicrocomputer. Design cbjectives are to 
provide an accurate navigation aid for general aviation in 
the $1000 price range. 


riTRCDUCTIOn 

Ongoing work at the Avionics Engineering Center, Ohio 
University, is aired at developing a simplified Oreca 
Navigation System (CNS) in the SICOC price range for use 
in general aviation aircraft. A low-cost Cr.oga senscr 
processor and unigue hardware phase-processing techniques 
have been developed and reported. See Burhans [1], 
Chamberlin [2,3] and Lilley [4]. 



Receiver Fronf-End Tviodule MiLfor.’ •npoler Inlerforc fAodule 


Figure 1. Sunnary Block Diagram - Kicrccoinputer-Based Omega Receiver. 
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In a related effort, a ioflware-based ONS using an In- 
eupcnslve mlcrocomputi r oiid miniiTui:! receiver bargware Is 
now being cxplori-d. toHptiler simulations of the proposed 
ONS are umleewiiy which use a data base 1 ive O.ioga phase 
data collcctid with th.: Ohio University Oi:;cqa sensor and 
microcomputer inttrfacc modules. Once the navigation 
algorithms have been fully developed in FCRTRAfj language, 
they win be translated into microprocessor cooe arrt loaded 
into the dedicated Of.S microcoii.puter, first as volatile 
progranialng, subject to change, and later as read-only- 
memory programs. 

This paper describes the ffi.cga sensor module, micro- 
computer Interface module, the establishment of the data 
base, and concludes with an ex-i:ple of the algorithms 
being developed. Operational results and low-cost con- 
siderations arc ciiphasirtd thnughout. 

The authors acknowledge iViSA/Lar.gley Research Center 
which has supported this v.ork under Grant I.'0R-36-009-01 7 
for the application of VLF techniques to general aviation 
air transportation systems. 

0.MEGA RECEIVING SYSTEM 



Figure 1 Illustrates In general form the receiving 
system to be described in tnis paper. To minirl^e cost 
and complexity of installation, tee system is cesigned to 
share the sense antenna of an aircraft autcmatic direction- 
finder (ADF) receiver. 

The preamplifier module, located very near tne antenna, 
provides Omega signal to coaxial caole at Ior Imoeeance. 
Simultaneously, the prear.pl if icr outputs o signal appro- 
priate for the sense antenna input of t^o ACF receiver. 
Onega signals arrive at the roceivc-r front-e-d r.pc’uie 
which (1) powers the preanpl if ler through the signal 
cable, (2) arplifle* the Onega si grai further arc filters 
to provide approxim.ately a 30 Hz car.dwidtn and (3) orcduces 
output pulses coincident with ure.'a signal zero-crossings. 
After processing in the Microcomputer interface r.cdule, 
the resulting Cnega pnase data is passed to a ciaital tape 
recorder (in the sim.ulation cor.figuraticn) or to the micro- 
computer as input to tre sensor a>'d navigaticn algcrith.r.s. 


Figure 3. Cmega/ADF Preamplifier Performance. 


RECEIVER FRONT-END MODULE 

The front-end module, described In detail by Surhans (11 
provides additional gain, indication of signal strength 
and digital pulse outputs at Omega signal zero-crossings 
ihe unit consists of three integrated-circui t chips plus’ 
two ceramic filters providing a bandwidth of appro-slmatel v 
.0 Hz at 10.2 KHz. See Figure 4 for the resoonse curve 
and Figure 5 for a summary of frcnt-erd m.odule so»cifi- 
catiens. 

The preamplifier and frcnt-er.d module have been tench- 
and f iph.t-tested as reported by Wright [6]. C.-ega phase 
-a .a taken from t.he 30 Hz bandwidth system has been ana- 
lyzed and reported by Zervos [7], 


PREA.MPLIFIER MODULE 


MICROCOMPUTER INTERFACE MODULE 


The preamplifier utilizes a field-effect transistor 
circuit providing a gain of 2C decibels at the 10.2 KHz 
frequency. The phase shift at tnis frequency is adjusted 
for zero degrees. At the ADF output, zero phase shift is 
seen across'the ADF band, ard t.he net gain is approxi- 
mately -6 decibels. See Figure 2 for details of pre- 
amplifier circuitry and figure 3 for performance 
curves [5]. 





The microco.Tputer Interface nodule performs the func- 
tions of pnase detection and dat.a sampling, it generates 
microcomputer interrupts for reading the pnase data ir.*o 
the .■nicrocemputer. 



Figure 2. ftrega/ADF Preamplifier Circuitry. 


Figure 4. 10.2 KHz Filter-Amplifier Response. 
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Figure 5. Receiver Front-End Specifications Surjr.ary. 


As a special case of a nul tiple-frequency ONS, a func- 
tional description will be given for a 10.2 KHz single- 
frequency system. However, the described tecr.niques are 
directly extendible to the rul tlple-cnennel case by simply 
deriving one additional clock treouency per channel from 
the existing ccrr.tion refererce oscillator (TCXO). A further 
specialization erplcyed here Is the use of a 2° x 10.2 Khz 
■ 2.6112 KHz reference clock to obtain 8-bit data word 
lengths. Longer or snorter ».ord leroths can be obtained 
by using higher or lo./er reference clock frequencies. 

For design generality, the clock frequency nctdtc for 
the phase detector's reference ccunter is presently derived 
from a 5 KHz TCXO and a rate rultipr.er countdown cnain. 

A multiple-frequency system cculd utilize a cerron TCXO 
and multiple synthesizer chains. To obtain 2.6112 .KHz. the 
5 HHz Is multiplied by 0.52224 in five decades of binary- 
coded-declral rate multipliers. The rate multiplier pro- 
vides for any multiplication factor between 0.1 ard 0.9 
using BCD Inputs to set the output rate. The output con- 
tains the required nu.'ster of pulses in unit time, tut due 


The phase measurenent of each Omega zero-crossing with 
resoect to a local reference clock is performed in an open 
loop fashion as shown in Figure 6a. The functional des- 
cription of the phase detection process is as follows: A 

clock signal of frequency 2^ x 10.2 KHz * 2.6112 KHz is 
applied to an eight-stage binary ripple counter. The out- 
put of the first stage (LSB) of this binary divider is. 
therefore, a square wave of frequency 2' x 10.2 KHz. Simi- 
larly, the last stage output (USD) is changing at a 2^ x 
10.2 KHz rate. The eight counter outputs are used as data 
inputs to an 8-bit latch. 

The latch is cen^posed of eight edge-clocked D flip-flops 
which are simultaneously clocked by the inceming Cnega 
zero-crossings supplied from the receiver front-end. Since 
data is transferred from the 0 Inputs to the 0 outputs of 
the latch only when it is clocked by an Omega zero-cross- 
ing, an 8-bit word representing the phase of the incoming 
Onega edge with respect to the local clock appears at the 
latch output approximately 10.2C0 tim.es per second. 

Note that no attempt is made to make absolute phase 
measurements; only a measurement-to-neasurement relative 
phase Is obtained at this point. For exam.ole, if the 
C.-ega signal was noise-free (a stable 10.2 KHz signal), 
the latch would be decked at regular intervals of T • 
1/10.2 KHz • 9.8 X 10‘= seconds. This wculd result in the 
output remaining constant for each measurement (see Figure 
,’a). This constant number would represent the initial 
prase offset between the inceming 10.2 KHz signal and the 
i" X 1C. 2 K.HZ reference clock. Realistically, however, 
since Omega signals off-the-air are contaminated oy noise 
(prase jitter), tre Cnega zero-crossings used to clock tne 
data latch will be irregularly spaced in tii e. T.his re- 
sults in rsneonstant prase measurement numbers at the latch 
output, the short term variations in this prase measure- 
ment being caused by the noisy Otega signal (Figure 7b). 



Figure 6. Microcomputer Interface Module. 


originacpag^ 

OF POOR QUALH^ 


- 3 - 








Figure 7. Phase Detection Waveforris. 


DATA SAMPLING AND INTERRUPT GEflERATION 

There are several considerations which suggest sampling 
the phase measurements at a rate slower than the 1C. 2 KHz 
carrier rate. 

(1) The high degree of correlation of the Or.ega signal 
makes It unnecessary to use the phase measurement from 
every cycle of the carrier. 

(2) Sampling of the data provides, in effect, an amount 
of needed integration if the sampling rate is slow 
enough to undersanple short-term phase variations due 

to noise. 

(3) Using the data at the maximum available raw-data 
rate implies a prohibitively high interrupt rate for 
most riicrcccfflputers. 

The choice of a particular sampling rate is based on 
these considerations. 

In this particular hardware implementation, a sampling 
frequency of 100 Hz has been chosen for the following 
reasons: 

(1) A front-end receiver bandwidth of 30 Hz Implies 
that samples must be taken at greater than a 60 Hz 
rate for the sampled signal to contain all the informa- 
tion of the raw input signal (Nyguist Rate). 

(2) Previous experience with hardware processing 
schemes indicates that atout 10 bits (1024 saniples) 

of Integration is optimum for this receiver configura- 
tion [3J. A 100 Hz sampling rate accepts every 102'’’^ 
phase measurement, in effect, resulting in nearly 7 bits 
of Integration. 

(3) 100 Hz is an acceptable microcomputer interrupt 
rate, and it is easily generated iron the existing 
receiver timing circuitry. 

Figure 6b shows the derivation of the ICO Hz sampling 
frequency from cxistinn timing signals. The ripple 
counter's i!SB which is toggling .it a 10.2 KHz rate is 
further divided by 102 in a prrgrai:.'’able divider chain. 
Similarly, for a mul tiple-freguoncy systmi. inis same 
100 Hz could be used to sample sii'^ultarcously tnc otner 
data channels as v;ell. The resulting 100 Hz sguare wave 
Is then fed to a positive-edge trigger to generate a 
m1croco»:puter interrupt every 10 msec. 

To Insure that the data at the output of the latch 
1 $ stable when the microcomputer interrupt occurs 


(I.e,, to Insure that the pseudo- random Omega zero-cross- 
ings do not clock new phase measu 'ements through the 
latch while the data Is being read), the Omega zero- 
crossings are gated with the interrupt pulse as shown In 
the microcomputer Interface module block diagram in 
Figure 6. 

In keeping with the design objectives of low-cost, 
simplicity and small size, this single-frequency version 
of the microcomputer Interface module has been orass- 
boarded using ten TTL integrated circuit chips. Further 
plans call for a CKOS version of this module to be im- 
plemented using just four chips. 

DATA BASE 

Before attempting to use a ccnmercially available micro- 
computer dedicated to the ONE, it is desirable to simu- 
late the navigation system functions and develop efficient 
(and simplified) algorithms using a high level programming 
language and a multipurpose ccmputcr. The classical 
approach to system simulations is to use computer number 
generators to approximate the signal and noise inputs 
from the known statistics of both. This approach would 
Involve passing an order Harkov signal (highly cor- 
related Cmega) and V.’hite Gaussian noise plus irpulse noise 
(VLF atmospheric noise) through an order bandpass 
filter (narrowband receiver front-end). 

However, a much more powerful approach has been taken 
to Insure that the data for the navigation processor 
simulations is exactly the same cata that will ulti- 
mately be seen by the dedicated microcomputer. Since 
the Cmega sensor end microcomputer interface mcdules 
of the proposed O'iS have been fully developed in ha. .lare 
and flight tested, it has been possible to use these same 
modules that will supply the navigation processor with 
Onega data to collect actual off-air data for use in the 
navigation processor simulations. For the collection of 
ground and flight data the Omega sensor and microcom.puter 
interface modules have teen used to supply Omega pht-e 
data to an incremental magnetic tape recorder. The 
microcomputer Interrupt pulse is used as the tape recuixer 
"v.rite" pulse. 

Twenty-feur continucus hours of ground data and five 
hours of flight data have teen collected as a data base 
fer si.-.ulations.. The Cmega sensor and microcomputer 
interface modules and Kennedy tape unit were installed 
in an Ohio University DC-3 flight test aircraft and used 
fer airborne data collection. 

NAVIGATIO:i PROCESSOR SI.'iULATION 

One goal of this project is to Investigate the concept 
of a microcomputer-based ONS; that is, to develop a 
ccm.pletely software-based system recuiring as litt’e 
hardware as possible. Previously reported research [4] 
has resulted in a hardware sensor processor ircludirg a 
hardware Cinega synchronization scheme and correl,' ticn 
detector. The proposed software-based system would per- 
form these functions as well as Skywave Correct on. 

Position Location, Coordinate Conversion, and fackoit 
Display Generation via software routines. Following the 
ccnpletion of this study, a direct ccoTarison of hardware 
versus software performance will be possible A quantita- 
tive evaluation of the hardware/sof tv/are enrineerirg 
tradeoffs in terms of performance, cost, anJ physical size 
can be performed. 

An example of the OMS algorithms being developed is 
the automatic synchronization to the Cmog.i transmission 
pattern which is described below. 

AUTO.KATIC SYNCHRONIZATION 

The Omega system consists of eight stations transmitt- 
ing on three frequencies in a fixed, tine-ir.ul liplexed 
fom-at. A unique p.itlcrn is formed by the scheduled 
lengtn of each station's transmission which repeats 
every ten seconds (see Figure R). Thus, to use the 
pK.ise measurfments supplied by the microcomputer inter- 
face n;oduIe, some means n.ust be provided for synchronizing 
the ONS timing with the Omega tranwiission pattern. Only 
after the O'iS is "in sync" can the ph.ise of each station's 
signal be followed by tracking loops to give position- 
fixing information. 
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Figure 8. Oraega Transmission Pattern. 


It should be pointed out that precision sync (within 
100 msec) is unnecessary. A .2 second gap between each 
Iransmission burst allows for propagation delays from 
transmitter to receiver without signal overlap into 
another station's time slot. Therefore, sync within 
♦ ICO msec insures that only one station will be received 
during each of the eight time slots. 

A description of the automatic synchronization procedure 
developed and simulated at Ohio University follows, and 
a logic flow diagram is given in Figure 9. Although the 
routine is described here in terns of a sinqie-freguency 
system, the technique is applicable to a multiple fre- 
quency system as well. Performing the routine on multiple 
channels simultaneously would give the navigation pro- 
cessor a redundancy check. By averaging the results 
obtained for all channels, a best estimate of the sync 
point could be obtained. 
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Figure 9. Automatic Svnc'r m za tion Flow Chart. 
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To establish sync we must obtain a transmission pattern 
from the live data and shift the entire ten second frame 
until the live pattern coincides with the unique Onega 
transmission pattern stored in the computer. The live 
transmission pattern can be established by determining 
the amount of correlation between consecutive lO-msec 
phase samples. High correlation indicates the presence 
of a strong signal, while the lack of phase coherence is 
indicative of a weak signal cr noise. The absolute value 
of phase measurements is meaningless during the sync pro- 
cess; only the phase correlation from sample to sample is 
useful in determining the live transmission pattern (i.e. 
in determining whether the sample point is signal or 
noise). 

An Omega phase measurement is supplied by the micro- 
computer interface module every 10 msec (100 Hz sampl- 
ing rate). Starting at a random point in tine and tak- 
ing data for ten seconds (one complete Omega frame) 
yields 1000 phase samples. This is all the information 
needed to accomplish sync: One ten second record con- 

taining 1000 bytes of phase data, each byte being 8-b1ts 
long. To simulate a random ONS start time, a call to 
a subroutine returns two random numbers. These represent 
the number of records and bytes in the data base tape to 
be skipped to get to the random starting point. Begin- 
ning at this point on the data base tape, 1000 phase bytes 
are read. 

The Omega signal is contaminated with noise, causing 
a series of measurements to result in a dispersion of 
phase points rather than a constant value. If the phase 
of a signal with respect to the receiver clock is near 
the 0 to 255 measurement extrem.es (near the ’’edge" of 
a cycle), the measured value may "bobble" between the 0- 
255 extremes. To rem.ove this bobble each data point is 
"mirrored" about t*'e 50” full scale value (mirrored about 
127). Points at X ' 255 are thus reflected to 255 - x = 0. 
Since we are interested only in the phase correlation 
and not the absolute value of phase, the mirroring pro- 
cess merely puts the data into a bettor form from which 
the correlation can be determined. 

A simple lOutine yielding a ...easure of the point-to- 
point correlation is to take the difference cf eech two 
successive phase points. That is, subtracting each 
measurement from the previous yields ICOO difference 
values: each difference is inversely proportional to 

the correlation between two sample slots. 

By comparing each difference value to a predeter- 
mined "threshold constant", a decision can be made as to 
whether signal or noise exists in each 10-msec sample slot. 
However, it may be desirable to smooth or average the 
difference values over several sample slots before com- 
parison to the threshold constant is made. If the differ- 
ence is less than the threshold, a "1" is stored in the 
sample slot (low difference — high correlation — > signal 
present). If the difference is greater than the threshold, 
a "0" is stored in the slot. After all lOCO comparisons 
have been performed, a pattern of 1000 I's and O's re- 
presents the signal and noise sample slots, respectively. 

Note that the 1000 bits can be stored in 125 B-bit 
bytes of memory (an amount of storage not prohibitive for 
simple microcomputer systems). Also bear in mind that 
virtually the entire microcomputer will be dedicated to 
the autosync routine during the sync process, since no 
Omega navigation information is available until after sync 
is complete. 

The 1000 sample bits are next compared to the stored 
Omega pattern in a bilevel correlation process. Use is 
made of both time slot and transmission gap information; 
whenever the sample slot binit is equal to the stored slot 
binit, the correlation counter is incremented by one. 

After all 1000 slots have been compared, the stored 
pattern is shifted by one slot and the comparison process 
is repeated. This shift and compare iteration is con- 
tinued until all loco possible stiiftcd patterns have been 
tested. The number of shifts necessary to obtain the 
highest bilevel correlation value represents the starting 
point for the "A" time slot. 

Having gone to the sync point in the data. Figure 10 
shows the resultant 1000 byte record plotted 10 points per 
line. The higher the signal-to-noise ratio of the station 
being received in each time slot, the less dispersion there 
is in the phase points plotted. 


For a confidence measure, the entire sync process can be 
performed several times and a best estimate of the sync 
point computed. If multiple frequencies are available, the 
process can be performed on each channel and the sync 
points averaged (or a sync point can be considered accept- 
able only if it results on 2 out ot 3 channels, etc.). 
Another possibility is to save the second highest correla- 
tion value as well as the highest, and consider the sync 
point acceptable only if the highest value is better than 
the second highest by a predetermined amount. 

Nothing has yet been said about operator assistance dur- 
ing the sync process. However, the operator's a priori 
knowledge of which stations are on the air and which ones 
will be well-received in his particular geographic area 
can be incorporated into the creation of the unique trans- 
mission pattern stored in the computer. This facilitates 
the bilevel correlation process, resulting in a more 
accurate sync point. 

Besides the sync point information made available by the 
euto-sync routine, a signal-to-noise ratio measurement 
can be derived very sinply from the bilevel sample slot 
oattern. After the sync point has been determined, the 
first 90 binits in each time slot can be summed; the sum 
obtained for each time slot is proportional to the SNR in 
each time slot. This information can be used by the 
navigation processor to automatically determine which 
stations are best for navigation (in the SNR sense). 

After the sync point has been determined and the ONS 
timing synchronized with the transmission pattern, the 
125 bytes of sample slot binits need not be saved unless 
it is desired to check the sync point at later tires to 
insure that sync is maintained. Several methods for 
periodically checking the sync point are now being in- 
vestigated. 

In sunrary, the technique presented requires very simple 
microcomputer instructions (compare, subtract, and shift); 
uses very little storage (125 bytes naxirum); usually 
accomplishes sync in less than 20 seconds (depending on 
the number of redundancy checks desired). The routine has 
demonstrated acceptable accuracy (within ♦ 100 msec) using 
the ftrega data base described earlier. Continued research 
utilizing decision theory techniques is ongoing to deter- 
mine the threshold constant in an optimal fashion. Several 
types of data smoothing routines have been used to obtain 
the smoothed difference values; a 9-point straight averag- 
ing filter seems to give best results. 

FUTURE ACTIVITY 

The automatic synchronization algorithm is only the 
beginning. Immediately following acceptance of the pre- 
liminary sync algorithm, work will begin on the Omega 
tracking loop software, to be followed by navigation 
algorithms and display drivers. Current plans call for 
physical implementation cf the microcomputer after simula- 
tions of the sensor processor functions are complete. 

Flight tests of the software sensor processor are planned 
for summer 1976. 
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Figure 10. One Record of Omega Phase Data (1000 Bytes) 
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